Theory of 2D transport in graphene for correlated disorder 
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We theoretically revisit graphene transport properties as a function of carrier density, taking 
into account possible correlations in the spatial distribution of the Coulomb impurity disorder in 
the environment. We find that the charged impurity correlations give rise to a density dependent 
graphene conductivity, which agrees well qualitatively with the existing experimental data. We 
also find, quite unexpectedly, that the conductivity could increase with increasing impurity density 
if there is sufficient inter-impurity correlation present in the system. In particular, the linearity 
(sublinearity) of graphene conductivity at lower (higher) gate voltage is naturally explained as 
arising solely from impurity correlation effects in the Coulomb disorder. 



One of the most studied properties of graphene is its 
electrical conductivity as a function of the applied gate 
voltage vifhich translates directly into the carrier density 
(n) dependent conductivity a{n) [ij. The functional de- 
pendence of <y{n) at low temperatures contains informa- 
tion [l[ about the nature of disorder in the graphene envi- 
ronment giving rise to the dominant resistive carrier scat- 
tering mechanism. Although there is a well-accepted the- 
ory [l| for graphene transport involving an interplay be- 
tween long-range charged impurity and short-range disor- 
der scattering, the theory is not universally accepted and 
cannot explain all experimental observations, indicating 
the possibility of important missing ingredients 

In this work, we provide a qualitatively new theory 
for the cr(n) properties of graphene and introduce a new 
physical explanation for the experimental observations, 
i.e., we explain why cr[n) ~ n for 'small' or 'interixiedi- 
ate' n and cr(n) ^ constant for 'large' n, with a smooth 
nonlinear crossover between the two asymptotic behav- 
iors. We also provide theoretical results for CFmin, the 
graphene minimum conductivity at the Dirac point, us- 
ing our new theory. We concentrate on the nature of the 
underlying static disorder limiting graphene transport in 
currently available samples where phonon scattering ef- 
fects are relatively weak (compared with disorder scat- 
tering) even at room temperatures Q. The quantitative 
weakness of the electron-phonon interaction in graphene 
gives particular impetus to a thorough understanding of 
the disorder mechanisms limiting graphene conductivity 
since this may enable substantial enhancement of room 
temperature graphcne-based device speed for technolog- 
ical applications as disorder remains the primary resis- 
tive mechanism limiting graphene transport even at room 
temperatures. Therefore, a complete understanding of 
the disorder mechanisms controlling cr{n) in graphene at 
T = is of utmost importance both from fundamental 
and technological prospcctives. 

The most important features of the experimentally ob- 
served cr(n) in graphene are: (1) a nonuniversal 
sample-dependent minimum conductivity a{n « 0) = 



CTmin at the charge neutrality point (CNP) where the av- 
erage carrier density vanishes; (2) a linearly increasing, 
a{n) (X n , conductivity with increasing carrier density on 
both sides of the CNP upto some sample dependent char- 
acteristic carrier density; (3) a sublinear a{n) for high 
carrier density, making it appear that the very high den- 
sity cr(n) may be saturating. 

A successful model [J, 0- 12 1 for diffusive graphene car- 
rier transport incorporates two distinct scattering mecha- 
nisms with individual resistivity pc and ps , arising respec- 
tively from the long-range Coulomb disorder due to ran- 
dom background charged impurities and static zero-range 
(often called "short-range") disorder. The net graphene 
conductivity is then given by a = p^^ ~ {pc + Ps)^^ ■ It 
is easy to show that 0, 0-113 Pc ^/n and ps ~ constant 
in graphene, leading to cr(n) going as tT(rt) = ri[[A + Cn) 
where the constants A and C are known [l[ as func- 
tions of disorder parameters; A, arising from Coulomb 
disorder, depends on the impurity density (rii) (and also 
on their locations in space) and the background dielec- 
tric constant (k) whereas the constant C, arising from 
the short-range disorder P, [ll| , depends on the strength 
of the white-noise disorder characterizing the zero-range 
scattering. The relation a{n) = n/{A + Cn) explains 
the observed cr(n) behavior of graphene for n ^ since 
a{n < A/C) n, and a{n > A/C) ^ 1/C with a{n) 
showing sublinear (C + A/n)^^ behavior for n ~ A/C. 

The above-discussed scenario for disorder-limited 
graphene conductivity, with both long-range and short- 
range disorder playing important qualitative roles at in- 
termediate {rii < n ^ ^/C) and high (n > A/C) carrier 
densities respectively, has been experimentally verified 
by several groups [5|-|8|. There is, however, one seri- 
ous issue with this reasonable scenario: although the 
physical mechanism underlying the long-raiige disorder 
scattering is experimentally established [l|, y, @| to be 
the presence of unintentional charged impurity centers 
in the graphene environment, the physical origin of the 
short-range disorder scattering is unclear and experimen- 
tally obscure. Point defects (e.g. vacancies) are rare 
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in graphene producing negligible short-range disorder. 
There have also been occasional puzzling conductivity 
measurements [e.g., Ref. [l^l reported in the literature 
which do not appear to be easily explicable using the 
standard model of independent dual scattering by long- 
and short-range disorder playing equivalent roles. 

In this Letter we propose an alternative physical 
model for understanding disorder-limited crin) behav- 
ior in graphene. The model is simpler (and therefore, 
more appealing) than the standard model of independent 
dual disorder mechanisms because it requires only the 
long-range Coulomb disorder associated with the back- 
ground charged impurities eliminating completely the 
ad hoc short-range disorder necessary for explaining the 
high-density nonlinearity in a{n). Our model, therefore, 
eliminates the undesirable feature of the standard model, 
namely, no adjustable short-range scattering term with 
unknown physical origin needs to be arbitrarily added to 
the problem in order to explain the observed high-density 
subline ar a{n). 

The key to our model is the inclusion of some spatial 
correlations in the distribution of the charged impurity 
locations in the system, i.e., the charged impurities are 
no longer considered to be completely random spatially. 
Some impurity correlations are perfectly reasonable to 
assume since much of the fabrication and processing of 
graphene is done at room temperature (and in fact, of- 
ten thermal and/or current annealing is used in sample 
preparation), which is expected to lead to actual diffu- 
sion of the impurities producing an annealed, at least 
partially, correlated impurity configuration rather than 
a quenched uncorrelated random one. We show that 
the single assumption of impurity correlations, defined 
through a correlation length scale parameter rg, is suf- 
ficient to explain the qualitative features of the experi- 
mental cr(n) behavior using only disorder scattering by 
background charged impurities. 

To calculate the impurity correlations we use Monte 
Carlo simulations carried out on a 200 x 200 triangular 
lattice with 10^ averaging runs, periodic boundary con- 
ditions, and a lattice constant Oq = 4.92A which is two 
times the graphene lattice constant since the most closely 
packed phase of impurity atoms (e.g. K as in Ref. @) on 
graphene is likely to be an m x m phase with m = 2 
for K IJ]. Correlations are automatically introduced by 
virtue of the random positioning of the impurities at lat- 
tice sites with the correlation length ro < ri ~ {nrii)"-^/^. 
Our correlation model is physically motivated with the 
reasonable underlying assumption that two impurities 
cannot be arbitrarily close to each other (as they can 
be in the unphysical continuum random impurity model, 
where rg =0), and there must be a minimum separation 
between them. A reasonable continuum approximation 
to this discrete lattice model is given by the following 
pair distribution function g{r) (r is a 2D vector in the 
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FIG. 1; (a) Density plot of structure factor S(q) obtained 
from Monte Carlo simulations for Ui = 0.95 x 10^^ cm~'^, ag — 
4.92 A and ro — 5ao. (b) Structure factor S'(q) using Eq. ^ 
(solid line) and Monte Carlo simulations. Dot-dashed and 
dashed lines show the Monte Carlo results for two different 
directions of q from a;-axis, 6^ = and 8 = 30°, respectively. 
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for the impurity density distribution. Even though 
Eq. ([1]) is only an approximation the basic idea of a 
length scale tq defining the spatial impurity correlations 
is physically sound (with tq = for the purely ran- 
dom case). Impurity correlation effects enter the trans- 
port theory through the structure factor S'(q), given by 
S{q) ~ 1 + Hi J cfre^'^''^[g{r) — 1]. For uncorrelated 
random impurity scattering, as in the standard theory, 
g{r) = 1 always, and S{q) = 1. With Eq. ((!]), we have 



Siq) 



2nni — Ji{qrn) 
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where Ji{x) is the Bessel function of the first kind. Fig.[T] 
(a) shows the structure factor 5(q) obtained from the 
Monte Carlo simulations. Fig. [IJb) shows >S'(q) for both 
the random Monte Carlo realistic numerical model and 
the simple continuum analytic approximation [Eq. ([2])]. 
It is obvious that the analytic approximation captures 
well the essential features of the full numerical Monte 
Carlo simulation. 

The graphene carrier conductivity due to scattering by 
screened Coulomb disorder can now be calculated taking 
into account the impurity correlations, leading to ct = 
{e'^/h){gEpT)/{2h), where Ep is the Fermi energy, g = 4 
is the total degeneracy of graphene, and the transport 
relaxation time r is given by, [15 
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where vp is graphene Fermi velocity, kp the Fermi 
wavevector (kp ~ Ep /{hvp)), and the graphene fine 
structure constant (ts = e'^/{hvpti)). For uncorrelated 
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FIG. 2: (color online) The carrier density for a single disorder 
realization obtained from the TFD theory (a) for the uncor- 
related case and (b) ro = 10 ao with Ui = 0.95 x 10^^ cm~^. 
Carrier probability distribution function P{n) are shown in 
(c), (d), (e) for (n) = 0, 1.78, 7.7 x 10^^ cm'^ respectively. 
In (f) the ratio n-rms/rii is shown as a function of ro/n for 
m = 0.95 X 10^^ cm~^, solid lines, and m = 4.8 x 10^^ cm~^, 
dashed lines. We use (n) = 7.7, 3.14, 0.94, x 10^^ cm"^ for 
the solid lines (from top to bottom) and (n) — 8.34, 4.10, 1.7, 



random impurity scattering, tq = 0, g(r) = 1, S'(q) = 1, 
we recover the standard formula for Boltzmann con- 
ductivity by screened random charged impurity centers 
(lH [l2| . In addition to scattering charge impurities in- 
duce strong carrier density inhomogeneities in graphene, 
especially close to the CNP, that must be taken into 
account in the transport theory. To characterize these 
inhomo gen eities we use the Thomas-Fermi-Dirac, TFD, 
theory |16j assuming that the impurities are placed in 
a 2D plane at a distance d = 1 nm from the graphene 
layer. Fig. [2] (a), (b) show the carrier density profile for a 
single disorder realization for the uncorrelated case and 
correlated case (rp = 10 ao) for rii ~ 0.95 x 10^^ cm^^. 
We can see that in the correlated case the amplitude 
of the density fluctuations is much smaller than in the 
uncorrelated case. The TFD approach is very efficient 
and allows the calculation of disorder averaged quanti- 
ties such as the density root mean square, rirms, and the 
density probability distribution P{n). Figures [2] (c) , (d), 
(e) show P(n) at the CNP, and away from the Dirac point 
{rii = 0.95 X 10^^ cm~^). In each figure both the results 




FIG. 3: Calculated cr{n) with S{q) obtained from the Monte 
Carlo simulations, symbols, and ^(q) given by Eq. El solid 
lines for (a) Ui = 0.95 x 10^^ cm"^ and (b) rii = 4.8 x 10^^ 
cm~^. (c) and (d) show the results for a{{n)) obtained from 
the EMT. The value used for n; in (c) and (d) is the same as 
in (a) and (b) respectively. The insets in (c) and (d) show the 
value of amin as a function of ro/n. In (e), the resistivity p is 
shown as a function of impurity density m for different carrier 
densities with ro = 5ao. (f) The relationship between n/ro 
and \/nro where the conductivity is minimum. The dashed 
line is obtained using Eq. [S] 



for the uncorrelated case and the one for correlated case 
are shown. P(n) for the correlated case is in general over- 
all narrower than P(n) for the correlated case resulting 
in smaller values of ?irms as shown in Fig. [5] (f) in which 
nrms/ni as a function of ro/r; is plotted for different val- 
ues of the average density, (n) , and two different values of 
the impurity density, ni = 0.95 x 10^^ cm^^ ("low impu- 
rity density" ) for the solid lines, and rii = 4.8 x 10^^ cm~^ 
("high impurity density") for the dashed lines. 

We now present our results for the conductivity. The 
integral in Eq. ([3]) can be calculated analytically for 
"small" kp by expanding S{x) in the integrand giving: 



where, A 



cr(n) = An [l — a + Ba^n/n 



(4) 
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G2(rs)/ (2G'i(rs)). Note a < 1 in our model. The di- 
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mensionless functions Gi,2{rs) arc given by, Gi{x) = 

% + 37ra;2 + 40a;3[l-7ra;+|(5a;2-l)5(x)], where g{x) = 
scch~^{2x)/\/T^^4x^ for X < i and sec~^ {2x) / \/ Ax'^ - 1 
for a; > 5. Eq. (|4]) indicates that for small n, a{n) ^ 
An{l — a)~^, and for large n, cr{n) (1 — ric/n) where 
ric ~ {1 — a)ni/ {Ba?) ^ 0{l/nir'^). The crossover den- 
sity nc, where the sublinearity (n > n,.) manifests it- 
self, increases strongly with decreasing tq. This gener- 
ally implies that the higher mobility annealed samples 
should manifest stronger nonlinearity in cr(n), since an- 
nealing leads to stronger impurity correlations (and hence 
larger rg). This is exactly the experimental observation. 
While the resistivity within the standard random model 
increases linearly in n^, Eq. (|4]) indicates that the resis- 
tivity could decrease with increasing impurity density if 
there is sufficient inter-impurity correlations present in 
the system. This is due to the fact that, for fixed tq, 
higher density of impurities are more correlated causing 
S'(q) to be more strongly suppressed at low q. This is 
easy to see in the case in which ro = ao and n,; so high 
that ri = ro- In this extreme case the charge impurity 
distribution would be very correlated, indeed perfectly 
periodic, and the resistance, neglecting other scattering 
sources, would be zero. For each value of ro and carrier 
density n, the maximum resistivity is found to be at 

r.Jro = y^2(l-7rBnr2). (5) 

In Figs. ISlJa) and (b), we show calculated a{n) using 
different values of the impurity correlation parameters 
(ro) and 5(q) given by Eq. ([2]) and Monte Carlo simu- 
lations. The comparison between the two results shows 
that the analytic continuum correlation model is quali- 
tatively and quantitatively reliable. It is clear that, for 
the same value o/ro, the dirtier (cleaner) system shows 
stronger nonlinearity (linearity) in a fixed density range 
consistent with experimental observation since the larger 
impurity density rii of the dirtier system allows, in prin- 
ciple, for stronger correlation effects to manifest itself 
due to the fact that the crossover density Uc is smaller 
for larger rii. To describe the transport properties close 
to the CNP and take into account the strong disorder- 
induced carrier density inhomogeneities we use the ef- 
fective medium theory (EMT) [1^. Fig. |3] (c) and (d) 
show the EMT results for a{n). The insets in Fig. |3] (c) 
and (d) show the dependence of amim on the size of the 
correlation length Tq. Umim increases slowly with tq for 
I'o/fi < 0.5, but quite rapidly for ro/r^ > 0.5. Finally, 
Fig. ^e) shows that the resistivity (1/cr) is highly non- 
linear as a function of impurity density and the optimal 
fi/ro at which the conductivity is minimum [Fig. [3](f))]. 

The results shown in Fig. [3] strikingly demonstrate 
the full power of the impurity correlation model as it 
clearly produces the observed experimental behavior with 



strong sublinear behavior for stronger impurity correla- 
tions (i.e. larger ro). Annealing leads to stronger cor- 
relations among the impurities since the impurities can 
move around to locate to equilibrium sites, thus enhanc- 
ing ro , which strongly suppress the crossover carrier den- 
sity nc{^ "1^0 thus increasing the overall nonlinearity 
of a{n). In addition the theory explains the observed 
strong nonlinear a{n) in suspended graphene (3] where 
the thermal/current annealing is used routinely. Finally, 
graphene on hexagonal BN is likely to have significant 
correlations in the impurity locations imposed by the 
similarity between graphene and BN lattice structure. 
This implies stronger nonlinearity in the cr{n) dependence 
for graphene/BN system as has recently been observed 
experimentally [l7| . Although we have used a minimal 
model for impurity correlations, using a single correlation 
length parameter rg , which captures the essential physics 
of correlated impurity scattering, it should be straight- 
forward to improve the model with more sophisticated 
correlation models if experimental information on impu- 
rity correlations becomes available 

In summary, we provide a novel physically moti- 
vated explanation for the observed nonlinear behavior 
of graphene conductivity by showing that the inclusion 
of spatial correlations among the charged impurity loca- 
tions leads to a significant sublinear density dependence 
in the conductivity in contrast to the strictly linear-in- 
density graphene conductivity for uncorrelated random 
charged impurity scattering. The great merit of our the- 
ory is that it eliminates the need for an ad hoc zero-range 
defect scattering mechanism which has always been used 
in the standard model of graphene transport in order 
to phenomenologically explain the high-density sublinear 
behavior. Even though the short range disorder is not 
needed to explain the sublinear behavior in our model 
we do not exclude the possibility of short range disorder 
scattering in real graphene samples, which would just add 
as another resistive channel with constant conductivity. 
We mention that a recent experimental work[l3| reports 
graphene transport data in remarkable agreement with 
the theory developed herein. 
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